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We study the QCD phase diagram at nonzero baryon and isospin chemical potentials using the 
l/N c expansion. We find that there are two phase transitions between the hadronic phase and the 
quark gluon plasma phase. We discuss the consequences of this result for the universality class of 
the critical endpoint at nonzero baryon and zero isospin chemical potential. 



I. INTRODUCTION 



In order to understand neutron stars, the early Universe, and heavy ion collision experiments, it is necessary to better 
grasp the physics of strong interacting matter in extreme conditions. Therefore the study of the QCD phase diagram 
at nonzero temperature and densities is very important. The nonperturbative lattice simulations that successfully 
addressed problems at zero densities can be used at nonzero isospin density but not at nonzero baryon density 
because of the so-called "sign problem". One has therefore to rely on novel approaches to study QCD at nonzero 
baryon density 0, H S H, which corresponds to the most important physical situations. In particular, these new 
methods have been used to study the critical temperature that separates the hadronic phase from the quark gluon 
plasma phase. However, these studies are valid only at small chemical potentials. 

In a previous article, we have successfully used the l/N c expansion of QCD to explain some key properties of the 
critical temperature that separates the hadronic phase from the quark gluon plasma phase 0. In the present work 
we shall extend our study to the general case \x u ^ fid and investigate the chiral phase transitions. Several models 
have shown that the phase diagram might be qualitatively altered in this case 000]. We shall also describe the 
consequences of our results for the universality class of the critical endpoint. We shall restrict ourselves to phase 
transitions between the hadronic phase and the quark gluon plasma phase, i.e. to situations with /i/ < m^, avoiding 
conditions where the ground state becomes a pion superfluid atT = o[l(l[ll0[l3. 

II. CRITICAL TEMPERATURES 

We use the usual l/N c expansion with 't Hooft's coupling [Til ]. As explained in detail in IEII3, the 1/JVc expansion 
of the pressure reads: 

1 Nf 1 
p(T, {m f , n f }) = N 2 C (p (T) + — X>( T > m /> <"/) + °(jva)) - W 

where / are the different light quark flavors. The sets of diagrams that contribute to the pressure are shown in Fig. 1. 
In a finite volume, the specific heat will peak at the transition between the hadronic phase and the quark gluon 




FIG. 1: Sets of diagrams that contribute to the pressure. The symbol / denote the quark flavor and runs over all flavors. The 
first diagram is O(N^), and the second is 0(N C ). 

plasma phase. This peak might diverge in the thermodynamic limit, depending on whether there is agenuine phase 
transition or only a crossover. The specific heat can be derived directly from the pressure leading to [6j 

C V = N 2 C (co(T) + — £ c i( T ' <"/) + °( Jp)) ■ ( 2 ) 
c /=i 
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Therefore, the critical temperature, which can be obtained by solving dCy / dT\T c — is given by [||: 

N f 

= +0(^5. ( 3 ) 

c /=1 

where the reduced temperature t c = (T c — Tq)/Tq, with To the critical temperature at zero chemical potentials. 

The chiral susceptibilities, Xf — d 2 p/dm 2 , peak at the transition or crossover where chiral symmetry is partially 
restored (chiral symmetry cannot be completely restored since we restrict ourselves to the case mf ^ 0). The 
diagrams that contribute to \f U P to next-to-next-to-leading order are shown in Fig. 2. The 1/N C expansion of the 




FIG. 2: Sets of diagrams that contribute to the chiral susceptibility Xf- The crosses are insertions of the mass operator, / 
is the specific quark flavor associated with x/, and g runs over all quark flavors. The first diagram is 0(N C ), the second and 
third are 0{N°). 

chiral susceptibilities can be expressed as 

1 Nf 1 

Xf = N c (xo(T,f4) + — (N fX i(T,$)+Y2(x2(t4)+X 3 (»f» g ))+xM))+°(w))> ( 4 ) 

C 9=1 c 

where xo comes from the first diagram in Fig. 2, Xi.2,3 come from the second diagram in Fig. 2 with the ^(-dependence 
taken respectively from the outside quark loop only (xi), the inside quark loop only (X2), and both the inside and 
outside quark loops (X3), and X4 comes from the third diagram in Fig. 2. In the equation above, we have used 
that p(fi u ,[id) — p{—p u . —fid) because of CP, and thus that Xfil^uiHd) = Xf(~^u, —^d), and that for equal masses, 
Xu(Hu,Hd) = Xd(^d, P-u)- If we restrict ourselves to m u = mj, notice that the that these latter two properties imply 
that Xu = Xd for [i u = ±Hd, but that Xu 7^ Xd hi general since \j? u ^ /j, 2 ,. Therefore, since the critical temperature for 
the restoration of chiral symmetry is defined by a peak in x/, each flavor might have a different critical temperature 
when /4 ^ l4 - 

The critical temperature that corresponds to the restoration of chiral symmetry, T/, can be obtained from 
dxf j 9T\t, — 0. Therefore, we find that the reduced temperature for chiral symmetry restoration, tf = (Tf — Tq)/Tq, 
is given by 

1 Nf 1 

tf = t (m/) + w (Nfni^ + J^^D+n^ffMg)) +r 4 (/4)) +0(^2). (5) 

Lattice simulations have shown that there is only one phase transition when /i u = ±/i^, at least at small chemical 
potentials for N c = 2 and 3 0, H H [3 • 

We assume that this property holds at large N c . Therefore the critical 
temperatures t c and tf for each flavor should coincide, at least at small chemical potential. However, it is clear from 
@ and JSJ that t c and tf do not have the same 1/N C expansion: the chemical potential enters at a different order in 
the 1/N C expansion, and, at leading order, t c depends equally on all chemical potentials, whereas tf depends only on 
one chemical potential. As was shown in |6(], the qualitative properties of t c are in agreement with the lattice results. 
We therefore impose that t c — t u — td for fi u = ±fid as a function of Nt and N c in the l/N c expansion. These 
constraints imply that 

T (p 2 ) = 

r 3 (p 2 ) = r z (-n 2 ) (6) 
T 4 {p 2 ) = 

f(p 2 ) = r 1 ( M 2 )+r 2 ( M 2 )+r 3 ( Ai 2 ). 

Notice that the same conclusions for QCD with an even number of flavors can be reached by using the dependence of 
t c and tf on the quark masses rather than on the chemical potentials. 
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In general, we expect that 71,3 7^ 0. We therefore conclude that critical temperatures for two different flavors will 
be related by 

N f 

tf-t g = —(NfinivPf) - n(v%)) + X^Mwm) - T 3 ((i g n h )j) + O(j^), (7) 

c h=l c 

and that the critical temperatures are sensitive to all chemical potentials at the same order in the 1/N C expansion. 
Therefore the restoration of the chiral symmetry takes place at different temperatures at /u„ 7^ fi d , i.e. at nonzero 
baryon and isospin chemical potentials. This is important since most experiments are precisely done in these condi- 
tions. A similar relation between the critical temperatures for different flavors can also been written as a function of 
the quark masses at zero chemical potentials. 



III. PHASE DIAGRAM 



We shall now analyze the consequences of this 1/N C analysis for the QCD phase diagram. We shall restrict ourselves 
to small positive chemical potentials, i.e. fi u , fid > 0. Possible generic phase diagrams resulting from the above analysis 
are sketched in Fig. 3. 
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FIG. 3: Generic phase diagrams at fixed T in the (fi u ,^d) plane in the 1/N C expansion. The chiral condensates are indicated 
where they are large. The solid curves are first order phase transitions, and the dotted curves are crossovers. The dots are 
critical endpoints, and the stars denote the points that become tetracritical endpoints at some temperature. See text for more 
details. 



In the phase diagrams presented in Fig. 3, the temperature is high enough so that the pion condensation phase is 
absent for fi u , fi^ > (or occupies at most a small portion of the phase diagrams near the fi u and fid axes that we 
will ignore here). The temperature in these diagrams is below the temperature of the critical endpoint at fi u = fid, 
T e [l7l ITiL Il9l]. Above this temperature, all the transitions become crossovers. The points denoted by stars in the 
phase diagrams of Fig. 3 become tetracritical endpoints at T — T e . The universality class of these points is that of 
the Z2 x Z2 Ising model, since there are two massless sigma modes that correspond to the divergence of both the 
baryon and isospin susceptibilities. 

It is possible that accidental circumstances might alter the generic phase diagrams presented in Fig. 2. If t\ = 
and T3 = 0, the chiral phase transitions would take place at the same temperature and chemical potentials for all 
flavors. Another accident could lead to the merging of the two critical endpoints in the second phase diagram of 
Fig. 3 into one single point at the corner of the region where all chiral condensates are large. In this case, the critical 
endpoint at T e would be in the universality class of the Z 2 Ising model, with only one massless sigma mode. However, 
we consider these latter two phase diagrams to be accidental in the l/N c expansion: They are accidents in the same 
way as To = is an accident. 

Similar phase diagrams as those presented in Fig. 3 can also be obtained from a slightly modified version of a 
Random Matrix model fA llflj. In the original model, at sufficiently high temperatures, the pion condensation phase 
never appears for fi u , fid > 0, IE 13 an d the flavor sectors decouple from each other. We use the Random Matrix 
effective potential and artificially add flavor- mixing terms to it: 

T-) 2 + (<7d - m ) 2 



RAIT 



(8) 



1 



— log 



(((al - {fil + an*) + T ) + 4TX)((^ - iti + a£) + T 2 ) 2 + AT 2 fi 2 ) + b{fi 2 u - fi 2 ) 2 a 2 u a 2 ) 



where a u ^d represent the chiral condensates. Depending on the choice of the artificially introduced parameters a and 
6, one can obtain either of the phase diagrams presented in Fig. 3. 
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IV. CONCLUSIONS 

In this article we have used the 1/N C expansion to study the QCD phase diagram at nonzero baryon and isospin 
chemical potentials, which corresponds to the most common experimental conditions. We have limited our study 
to the phase transitions from the hadronic phase to the quark-gluon plasma phase. We have also used constraints 
derived from lattice simulations, in particular that the chiral and deconfinement phase transitions take place at the 
same temperature for small baryon and isospin chemical potentials. For QCD with two flavors, we have found that 
there are two chiral phase transitions in general. The QCD phase diagram at nonzero baryon chemical potential is thus 
qualitatively different at zero and nonzero isospin chemical potential, as other models have predicted 0, S @ • This 
is important since most experiments are done at nonzero baryon and isospin chemical potentials. As a consequence, 
the universality class of the critical endpoint at nonzero baryon and zero isospin chemical potentials should be that 
of the Z2 x Z2 Ising model. The nature of the critical endpoint has observable consequences for heavy ion collision 
experiments |2l|. We conjecture that the real phase diagram resembles the second one in Fig. 3, because, following 
the analysis of Pisarski and Wilczek |2C|, the restoration of chiral symmetry for one flavor should be a crossover. The 
question of the deconfinement phase transition remains open. However, we conjecture that it will take place where 
the first chiral transition happens, since it is widely believed that the deconfinement phase transition cannot happen 
after a chiral phase transition. 

Finally, we want to comment on the possibility to test the predictions made above. First the properties of the sets 
of diagrams that enter into the chiral susceptibility and in particular the relations JSJ can be tested on the large- N c 
lattice approach outlined in • We stress that these tests can also be carried at zero chemical potential by studying 
these diagrams as a function of the quark masses. Second, the current lattice techniques used for fi u = fi^ can easily 
be extended to the more general case fj, u ^ fid- It might however be difficult to see the two chiral phase transitions on 
the lattice in this case, since these studies arc limited to small chemical potentials and that the separation between the 
critical temperatures should be proportional to hbijli/N c . Third, the 1/N C expansion also leads to different critical 
temperatures for the different flavors if they have different masses. Recent lattice results at zero chemical potentials 
have shown that the critical temperatures indeed differ for two light quarks and one heavier quark |23j , in agreement 
with the X/N c expansion analysis presented here. 
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